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SUMMARY 

An  analysis  of  the  temperature  distributions  in  a fluid  over  an 
oscillating  surface  with  heat  transfer  is  made  and  associated  heat- 
transfer  parameters  are  compared  with  those  for  the  case  of  conduction 
at  a stationary  surface  with  the  same  initial  temperature  potential.. 

It  is  found  that  the  heat  transfer  for  the  oscillating  surface  can  be 
considerably  different  from  that  for  conduction  alone.  The  effect  of 
the  surface  oscillations  on  the  thermal  state  of  the  fluid  is  studied 
by  means  of  average  static-  or  total -temperature  defects,  and  it  is 
demonstrated  that  the  oscillations  could  alter  the  fluid  temperature 
appreciably. 


INTRODUCTION 

The  increased  improvement  of  present-day  propulsion  systems  and 
the  development  of  new  propulsion  systems  have  posed  numerous  new  prob- 
lems in  the  field  of  heat  transfer.  Elucidation  of  unusually  high  heat- 
transfer  coefficients  which  are  apparently  encountered  in  unsteady  flows 
and  means  of  increasing  heat-transfer  coefficients  under  given  condi- 
tions are  greatly  desired.  As  a preliminary  attempt  to  gain  insight 
into  such  problems,  it  seems  worth  while  to  consider  the  heat-transfer 
aspects  of  the  classical  problem  wherein  the  fluid  motion  is  induced  by 
oscillating  a conducting  surface  axially  in  viscous  fluid.  In  particu- 
lar, the  effect  of  disturbing  the  equilibrium  (steady  state)  conditions 
after  the  periodic  motion  of  the  fluid  has  been  established  will  be 
studied.  The  temperature  distributions  in  the  fluid  are  determined  as 
exact  closed-form  solutions  of  the  energy  equation  pertinent  to  the 
problem  and,  hence,  related  heat-transfer  parameters  can  be  compared 
with  those  for  a stationary  surface  to  demonstrate  the  effect  of  the 
surface  oscillations.  Other  exact  solutions  of  the  energy  equation  for 
somewhat  analogous  problems  are  presented  in  references  1 to  3. 
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ANALYSIS 

The  laminar  motion  of  an  incompressible  viscous  fluid  induced  by 
the  axial  (longitudinal)  oscillations  of  a bounding  plane  surface  doubly 
infinite  in  extent  is  well  known  and  is  reported  in  numerous  sources 
(see  Schlichting,  ref.  4,  for  example).  The  equations  of  motion  for 
this  case,  assuming  constant  physical  properties,  reduce  to  the  class- 
ical heat -conduct ion  equation 


chi 


d2u 

Sy2 


(i) 


where  u is  the  velocity  component  parallel  to  the  surface,  t denotes 
the  time,  v is  the  kinematic  viscosity,  and  y is  the  coordinate  nor- 
mal to  the  surface.  (All  symbols  are  defined  in  the  appendix.)  The 
associated  boundary  conditions  are 


and 


u(0,t)  «•  U cos  nt 


u(-,t)  * 0 


(2) 


(3) 


where  U denotes  the  amplitude  and  n the  frequency  of  the  oscilla- 
tions. The  velocity  distribution,  after  sufficient  time  has  elapsed  for 
the  periodic  motion  to  be  established,  is  given  by 


u(y,t)  ««  U exp  (-  y^  cos  ^nt 

The  appropriate  energy  equation  is 

dT  d2! 

=*  a — ? 
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(4) 


(5) 


where  T denotes  static  temperature,  a is  the  thermal  diffusivity,  and 
Cp  is  the  specific  heat  at  constant  pressure.  Note  that  the  last  term 

on  the  right  of  equation  (5)  is  due  to  aerodynamic  heating  and  describes 
the  increase  in  fluid  temperature  caused  by  the  oscillations.  In  order 
for  the  problem  to  be  physically  reasonable,  however,  after  a long  per- 
iod of  time  the  system  should  be  in  a state  of  thermal  equilibrium  ^that 
is,  in  essentially  a steady  state)  and  it  is  thus  necessary  to  cool  the 
surface.  The  thermal  equilibrium  state  can  be  described  by  first  solv- 
ing equation  (5)  subject  to  the  boundary  conditions 


T(0,t)  =*  T. 
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(6) 
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and 

T(-,t)  - Tm  (?) 

that  is,  the  oscillating  surface  is  maintained  at  some  uniform  tempera- 
ture Tw.  It  is  unnecessary  here  to  specify  any  initial  conditions,  for 
they  will  not  influence  the  results  after  a sufficiently  long  time. 

After  the  explicit  form  of  the  nonhomogeneous  or  aerodynamic  heat- 
ing term  in  equation  (5)  is  determined  by  means  of  equation  (4),  a par- 
ticular solution  of  equation  (5)  was  found  to  be 

tp  - - ^ -» (-  VBgr  y)  [i  + e-(-r;-Pf^^]  (e) 

where  Pr  * v/a  is  the  Prandtl  number.  The  complementary  function  to 
be  added  to  equation  (8)  in  order  that  the  complete  solution  satisfies 
the  boundary  conditions  given  by  equations  (6)  and  (7)  is 


- erf 


(i^) 


(9) 


where  erf ( ) is  the  error  function  and  is  defined  as 

,P 


'f(p)  S ^ f exp(-r2)  dr 


so  that 

erf (0)  - 0 
and 

erf(-)  - 1 

The  complementary  function  is  synthesized  from  well-known  solutions  of 
the  heat  equation;  for  example,  the  term  containing  the  product  of  the 
exponential  and  trigonometric  functions  corresponds  to  the  solution  of 
the  heat -conduct ion  equation  for  the  case  where  the  temperature  on  a 
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stationary  surface  is  specified  to  be  a harmonic  function  of  time  (see 
ref.  5,  p.  47),  and  the  two  error-function  terms  combined  correspond  to 
the  solution  of  the  problem  of  unsteady  heat  conduction  to  or  from  a 
stationary  surface  due  to  a temperature  difference  or  potential  of 
u2pr  - 

T„  - 43 — - Tw>  Note  that  the  error-function  terms  describe  the  tran- 
5 

sient  heat  transfer  which  will  vanish  for  large  time.  The  temperature 
distribution  in  the  equilibrium  state  (that  is,  after  a long  time)  should 
be  independent  of  time  except  for  periodic  terms  and  is  given  by  the  sum 
of  equations  (8)  and  (9)  with  the  error-function  terms  vanishing  as 


The  relation  for  thermal  equilibrium  between  the  surface  temperature  Ty 
and  the  ambient  temperature  T„  can  be  found  by  evaluating  equation 
(10)  as  y -*  « and  is  given  by 
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U2Pr 

4CP 


(11) 


Thus,  as  was  anticipated,  it  can  be  seen  from  equation  (ll)  that  for 
thermal  equilibrium  the  surface  must  be  cooled  to  compensate  for  the 
aerodynamic  heating,  and  the  temperature  distribution  above  a surface 
oscillating  for  a long  time  in  a fluid  at  temperature  T'  is  given  by 


Now  that  the  fully  developed  state  of  the  fluid  is  completely  de- 
scribed (by  eqs.  (4)  and  (l2)),  consideration  can  be  given  to  the  prob- 
lem at  hand:  namely,  the  effect  of  the  oscillations  on  the  heat  trans- 

fer if,  at  some  time  (say  t = 0)  after  the  periodic  motion  (steady 
state)  is  established,  the  surface  temperature  is  fixed  at  some  temper- 
ature Ty  different  from  that  given  by  equation  (ll).  The  problem  to  be 
solved  to  answer  this  question  is  mathematically  identical  to  the  one 
for  thermal  equilibrium  except  that  Tw  replaces  Tw  in  the  boundary 
conditions  and  that  the  solution  of  the  present  problem  must  match  that 
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given  by  equation  (l2)  for  t =»  0;  that  is,  the  proper  initial  condi- 
tion must  be  satisfied  by  the  temperature  distribution  in  order  to  dem- 
onstrate the  transient  effects  properly.  It  was  previously  shown  that 
the  sum  of  equations  (8)  and  (9; 


satisfied  the  differential  equation  (eq.  (5))  and  the  boundary  condi- 
tions (eqs.  (6)  and  (7)),  where  Tw  is  now  replaced  by  Tw.  Evalua- 
tion of  equations  (l2)  and  (l3)  at  t = 0 shows  that  they  are  identi- 
cal so  that  the  proper  initial  condition  is  satisfied  and  thus  equation 
(13)  represents  the  temperature  distribution  over  an  oscillating  sur- 
face which  is  maintained  at  a constant  temperature  Tw.  Because  Tw 
is  a surface  temperature  different  from  that  for  thermal  equilibrium, 
the  solution  given  by  equation  (l3)  is  of  greatest  physical  significance 
only  as  long  as  the  transient  terms  are  important.  This  of  course,  is 
the  problem  of  primary  interest  herein. 

Although  the  solution  given  by  equation  (13 ) was  developed  for  an 
incompressible  viscous  fluid  with  constant  property  values,  it  is  equally 
valid  in  the  case  of  a compressible  viscous  fluid  if  the  boundary-layer 
assumptions  are  made,  if  the  Prandtl  number  and  the  product  of  p and 
p are  taken  to  be  constants,  and  if  y is  replaced  by  q where 


1 “ 


dy 


This  follows  because  the  von  Mises  transformation  (see  ref.  6,  for 
example)  under  these  assumptions  reduces  the  compressible  boundary- 
layer  equations  to  the  forms  of  equations  (l)  and  (5) . 


Any  pertinent  heat-transfer  quantities  can  be  obtained  from  equa- 
tion (l3) . However,  since  the  time  dependence  itself  is  usually  not  of 
primary  practical  importance,  the  temperature  gradient  at  the  wall 
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(which  is  a measure  of  the  heat  transfer)  will  be  averaged  over  a period. 
This  average  gradient  is  given  by 


where  tQ  indicates  the  beginning  of  any  cycle.  For  the  equilibrium 
state,  an  analogous  equation  is  obtained  from  equation  (l2): 


(15) 


This  average  temperature  gradient  is  due  entirely  to  the  aerodynamic 
heating  and  is  also  essentially  (that  is,  except  for  the  multiplicative 
constant  k,  the  thermal  conductivity  coefficient)  equal  to  the  work 
done  per  cycle  in  oscillating  the  surface.  The  latter  result  can  be 
independently  verified  by  computing  the  work  done  per  cycle  from  the 
shear  stress  on  the  plate. 


RESULTS  AM)  DISCUSSICN 

Now  that  the  temperature  distribution  and  its  gradient  at  the  sur- 
face are  known,  the  effect  of  the  surface  oscillations  on  the  heat 
transfer  can  be  obtained  by  comparing  any  appropriate  parameters  with 
the  corresponding  ones  for  the  case  of  pure  conduction  to  or  from  a 
stationary  surface  subject  to  the  same  initial  temperature  difference 
(T.  - Tw).  The  specific  formulation  of  the  parameters  depends  on  the 

particular  configuration  considered.  With  this  in  mind,  the  parameters 
considered  herein  will  fell  into  two  categories:  total  parameters,  that 

is,  those  pertaining  to  configurations  in  which  the  oscillations  sire 
inherent;  and  net  parameters  for  configurations  in  which  the  work  done 
in  oscillating  the  surface  is  assessed  to  the  system.  For  example,  the 
average  total  heat-transfer  coefficient  as  given  essentially  by  equation 
(l4)  can  be  compared  with  that  for  conduction  to  or  from  a stationary 
surface  at  an  initial  temperature  difference  (T„  - Tw).  The  temperature 
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distribution  for  the  latter  case  is  given  by  the  sum  of  the  first  and 
: . last  terms  of  equation  (l3)  (see  ref.  5,  p.  34)  and  the  associated  aver- 

j age  temperature  gradient  is  given  by  merely  the  last  term  of  equation 

I (14) . Therefore,  the  ratio  of  these  average  total  heat-transfer  coef- 
ficients is  given  by 


It  can  thus  be  seen  that  the  oscillations  could  lead  to  higher  average 
total  heat-transfer  coefficients  depending  on  the  sign  of  (T„  - Tw)  and 
on  the  relative  magnitude  of  the  two  terms  in  the  bracketed  part  of  equa- 
tion (16);  for  example,  in  the  limiting  case  just  after  the  temperature 
Tw  is  imposed  on  the  surface,  that  is,  tQ  -*•  0,  equation  (l6)  reduces  to 

\0  - - (VS  - »)  &7) 

so  that  for  the  case  where  the  ambient  temperature  is  greater  than  the 
surface  temperature  (T„  >TW),  and  Pr<  it  , then  R10>lj  that  is,  the 
oscillations  increase  the  average  total  heat -transfer  coefficient.  Con- 
versely if  Tw  > T„,  R10  would  be  larger  than  unity  if  Pr  > jt2 . 

It  should  be  reiterated  that  the  parameters  for  comparison  depend 
upon  the  configuration  considered.  Therefore,  if  in  the  specific  appli- 
cation the  net  energy  obtained  from  this  system  (rather  than  the  total 
as  before)  is  of  interest  then  the  total  heat  transfer  must  be  properly 
assessed  to  take  account  of  the  work  done  in  oscillating  the  surface. 
Hence,  the  net  heat  transfer  for  this  case  (as  represented  by  the  dif- 
ference between  eqs.  (l4)  and  (l5))  will  be  compared  with  that  for  con- 
duction to  or  from  a stationary  surface.  The  ratio  Rg  of  the  average 
net  heat  transfer  for  the  oscillating  surface  to  that  for  a stationary 
surface  can  therefore  be  written  as 


It  can  be  seen  from  equation  (l8)  that  the  average  net  heat-transfer 
coefficient  for  the  oscillating  surface  would  be  larger  than  that  for 
the  stationary  surface  (that  is,  Rg>l)  only  when  TW>T„,  and  it  would 
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be  lower  when  Tw<  T„.  For  a given  temperature  difference,  therefore, 
the  ratio  of  the  average  net  heat-transfer  coefficients,  as  given  by- 
equation  (18),  can  be  altered  by  changing  the  amplitude  of  the 
oscillations . 

If  the  primary  interest  is  not  in  the  net  heat  transferred  across 
the  surface  (as  it  would  be,  for  example,  for  regenerators)  but  is  in 
the  thermal  state  of  the  fluid  itself,  the  effect  of  the  oscillations 
can  be  studied  by  comparing  the  average  static  temperature  defect  de- 
fined as 


tn+  — 
I u n 


(T„  - T)dy  dt 


for  the  oscillating  surface  to  that  for  the  stationary  surface.  The 
ratio  of  these  defects  is  given  by 


R3=  W 


ifipr  .[pint  / 

" 4cp  (T„-Tw)  + 8ncp(T.-Tw)  ^ 2n  |_V 


0 + n) 


For  the  limiting  case,  tQ  -*■  0,  equation  (19)  reduces  to  ^ 

R » 1 7-^ r (l  - - iJPr)  (21) 

30  4Cp(T-  - Tw)  \ 8 y J 

It  can  thus  be  seen  that  unless  Pr>64/9  the  oscillations  will  not  be 
advantageous  for  lowering  the  fluid  temperature  for  the  case  when 
T >T„. 

A more  meaningful  parameter  to  examine  would  perhaps  be  the  ratio 
of  the  average  total-temperature  defect  with  the  oscillations  to  that 
for  conduction  to  the  stationary  surface.  The  average  total -temperature 
defect  is  defined  as  in  equation  (19)  except  that  T is  replaced  by  the 
total  temperature  Tt,  which  is  given  by 
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In  this  way  the  thermal  condition  of  the  fluid  independent  of  its  motion 
is  indicated.  This  ratio  is  given  by 


Thus,  the  ratio  R4  could  be  greater  than  unity  dependent  on  the  rela- 
tive orders  of  magnitude  of  the  terms  in  the  braces  in  equation  (22) 
and  in  the  limiting  case  as  tg  -►  0. 


R 


‘40 


1 - 


U2Pr 


4Cp(T. 


—)  1 + 
w' 


3 (l  - Pr) 

8 v** 


i 


(23) 


Thus,  in  this  case  if  T„ > T^  the  fluid  would  be  at  a lower  relative 
temperature  with  the  oscillations  (R40 > l)  only  if  Pr>9. 

An  analogous  discussion,  of  course,  follows  for  increasing  the 
fluid  temperature. 


CONCLUDING  REMARKS 

The  results  of  this  preliminary  analysis  of  the  temperature  distri- 
bution over  conducting  oscillating  surfaces  have  shown  that  the  heat 
transfer  associated  with  the  oscillating  surface  can  be  significantly 
different  from  that  for  conduction  between  the  same  initial  temperature 
difference  and  a stationary  surface.  In  addition  to  the  heat  transfer, 
the  thermal  state  of  the  fluid  was  studied  by  means  of  average  static - 
or  total-temperature  defects  and  it  was  found  that  the  surface  oscil- 
lations can  alter  the  thermal  state  of  the  fluid.  In  each  case,  the 
range  of  parametric  values  for  which  the  surface  oscillations  are  bene- 
ficial were  indicated. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  February  19,  1954 
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APPENDIX  - SYMBOLS 

The  following  symbols  are  used  in  this  report: 
specific  heat  at  constant  pressure 

P 


R 


10 


30 


40 


R, 

V 

R 

R4 

R4 

T 

Tt 
T, 

% 


w 


— P exp  (-r^)  dr 

Vn  Jo 


thermal  conductivity  coefficient 
frequency  of  oscillations 
Prandtl  number 

ratio  of  average  total  heat -transfer  coefficients 

ratio  of  average  total  heat-transfer  coefficients  evaluated  as 
tQ  -0 

ratio  of  average  net  heat-transfer  coefficients 

ratio  of  average  static-temperature  defects 

ratio  of  average  static -temperature  defects  evaluated  as 
tQ  -0 

ratio  of  average  total-temperature  defects 

ratio  of  average  total-temperature  defects  evaluated  as 
tQ  -0 

static  temperature 

equilibrium  (or  steady-state)  static  temperature 
total  temperature 
surface  temperature 

surface  temperature  for  steady  state 
static  temperature  in  ambient  fluid 
temperature  defect 
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t time 

tQ  time  at  the  start  of  an  arbitrary  cycle 

U amplitude  of  oscillations 

u velocity  parallel  to  surface 

y coordinate  normal,  to  surface 

a thermal  diffusivity 

Tj  von  Mises  coordinate  normal  to  surface 

H absolute  viscosity 

v kinematic  viscosity 

p density 


Subscripts: 

c complementary  solution 

osc  oscillating  surface 

p particular  solution 

sta  stationary  surface 
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